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On Lax-Wendroff-type time integration in high order 
finite difference schemes 
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Abstract — In this short note we present a new approach of 
classical Lax-Wendroff-type time integration for constructing 
high order finite difference schemes for solving the 
one-dimensional nonlinear hyperbolic conservation laws. This 
approach is based on one property of initial equations. The main 
advantage of this new approach is simplicity, because we use the 
degrees of Jacoban instead of its higher derivatives. For the first 
time finite difference schemes with M={ 3, 5, 7, 9,11} temporal 
and spatial order of accuracy were constructed. Results of 
selected classical test cases confirm accuracy and non-linear 
stability of these schemes in the cases of the nonlinear 
hyperbolic conservation laws. 

Index Terms —Lax-Wendroff-type time integration, finite 
difference scheme, high order of accuracy. 


I. INTRODUCTION 


In this paper we describe new approach to using 
Lax-Wendroff-type time integration [1] in constructing the 
high order finite difference schemes (FDS). We describe this 
approach in application to solving the one-dimensional 
nonlinear hyperbolic conservation law 

8f , dg{f) _ 
dt dx 

Here / is the function, which has to be found from (1), 
g(f) is the flux function. 


- = 0 . 


( 1 ) 


At present time, there are two main approaches in time 
integration [2]. The first approach is using ordinary 
differential equation (ODE) solver (usually this is 
Runge-Kutta solver) to solve system of ODE obtained after 
spatial discretization of (1). Second one is Lax-Wendroff 
approach, which is based on Taylor expansion on time. 

The first approach (Runge-Kutta-type time integration 
scheme) is used widely in theory of numerical methods for 
solving of partial derivative equations (PDE) and systems of 
PDEs [3]. The main advantages of this approach are its clear 
concept, simplicity in coding, and good stability property 
(especially if total variation diminishing (TVD) 
modifications of Runge-Kutta methods [r4] are used). The 
main lack of this approach is limit of the order of accuracy in 
TVD Runge-Kutta methods - they cannot exceed fourth order 
of temporal accuracy [5]. 

The second approach (Lax-Wendroff-type time integration 
schemes) is based on classical Lax-Wendroff FDS [1]. The 
main idea of it is in joint approximation of temporal and 
spatial derivatives. For this purpose, temporal derivatives in 
modified equation of FDS are replaced by spatial derivatives 
using initial PDE and its differential extensions. It leads to 
significant decreasing of the number of points in the 
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calculation stencil. Another good property is absence of 
limitations on the order of accuracy - virtually one can 
construct FDS of any desired order of accuracy. The main 
drawback of Lax-Wendroff-type time integration schemes is 
in their complexity, especially for systems of PDEs. In this 
case one has to operate with multidimensional tensors like 
8 n g(f)/df n [2]. Therefore at present time there are no 
Lax-Wendroff-type time integration schemes for equation (1) 
with higher than fourth order of temporal accuracy. 

The aim of this work is to present a way of constructing 
finite difference schemes for PDE (1) with 
Lax-Wendroff-type time integration with (any) high order of 
accuracy. The approach suggested below is deprived of 
complexity of existing Lax-Wendroff-type time integration 
schemes because it does not operate with the derivatives of 
the Jacobian of the PDE (1). At first, one property of the 
equation (1) is proved. After that, new Lax-Wendroff-type 
time integration scheme is described. Set of high order (up to 
eleventh order of temporal and spatial accuracy ) FDS for 
problem (1) is constructed. Finally numerical results show 
properties of constructed FDSs. 


II. P re liminary 

Together with equation (1) we consider its 
non-conservative form 

—+ A—=0. (2) 

dt dx 

Qg 

Here A(f) = — is the Jacobian of equation (1). We 

df 

prove one property of equation (1). 


Theorem. If f is solution of one-dimensional nonlinear 
hyperbolic conservation law (1) and all necessary derivatives 
of f and g(f ) exist, then following relations between 
temporal and spatial derivatives of f 



are correct. 

OWe conduct the proof to two steps. 

1. We prove first equality of set (3): 

d^d_f A 8 A 

dt 2 dx v dx) 

Its correctness follows from the sequence of equalities, 
presented below: 
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dt dx\ = 

dg = Adf \ 




d 

(*1 

dt{ 

dx) 

dx 

UJ 

--I 

A— 


A— 


v dt. 

J dx{ 

. dx. 


2. We prove, that if equality (3) is correct for k=n 

d n f r-'f, ^ 0 - 1)30 

dt n Sx"" 1 ^ ’ dx)' 

then it will be correct for k = n +1 

dt n+l 5x"l l ’ dx/ 

First, we examine expression for 5” +1 / /5t” +1 : 

_ L fel -L{ -^((-aY - 1 %]) - 

8t n+1 dt{dt n ) 4 ; dx)) 

■■SIM *<-*>-& 


Further, we consider the summands in the brackets in right 
part of (4): 


&g_ 

=-[ 


l=-l 

A—1 = - 


A—). 

dtdx 

dt ! 

,dt j 

1 5x1 

V St] 

dx\ 

V dx) 


Non-conservative form (2) of (1) was used intensively at 
chains (5, 6). 

To complete the proof of step 2, we substitute relations (5, 
6) into (4): 


d n+1 f 

dt n+l 



Finally, the correctness of theorem follows immediately from 
step 1 and step 2 on an induction. 0 


III. Construction of the Schemes 


In domain {0 < x < 1,0 < t < 7\} we consider uniform 
grid 

co xh = \ x i - ih I i = 1 ,N x ;tj -tn\n- 1, N t }, 
where X,h are uniform (for simplicity) temporal and 
spatial steps respectively. We denote the point value of any 
grid function cp(f n ,x,-) = cp {nx,Ui) as cp” . 

Following classical Lax-Wendroff-type time integration 
scheme [1, 2], we obtain semi-discrete form of initial 
equation (1). For this purpose we consider following 
expression for derivative df /dt: 


df_ fi n+1 ~fi n y t* d k+1 f 
dt ~ x f~l(k + l)\ dt k+1 


(7) 


Such representation follows from Taylor expansion of the 
finite difference (/) n+1 — around the point (f n ,x ( ). 


Substituting (7) in (1), we obtain semi-discrete form of initial 
equation (1): 

V' -/" ,8s f t‘ d^f 

t dx £f(k + l)! ar i+1 
Further, we replace temporal derivatives in (8), using (3): 


( 8 ) 


fi-f" 


dx ^+1)! dx k v ’ dx) 


(9) 


We remark, that equation (9) is completely equivalent to 
the initial one (1). The difference between them will arise, 
when finite number of terms in sum in (9) is taken into 
account. We have to consider (M -l) first terms in sum to 


achieve M order of temporal accuracy. 

We derive from equation (9) original two-step approach 
for constructing finite difference scheme with desired order 
of temporal and spatial accuracy M for solving (1): 

1. We reconstruct finite difference approximation of 
derivative dg/dx with M order of accuracy 


9; = Lm (#? ) 


Mill 

dx [ 


+ 0(/i M ). 


( 10 ) 


Hereinafter l} M is some difference operator, which 
approximates derivative d k /dx k with M spatial order of 


accuracy. 

2. We calculate on the new temporal layer using 
relation 

Af-t k +1 . , , 

Sr' = /"-W, - £ Jj—fj-pM-* (Hi) <4 <“> 

Any type of reconstruction (from simplest centered finite 
differences [6] up to modern used widely WENO 
reconstruction and its modifications [2,7,8]) may be used at 
first step for definition operator l} M . Using simplest centered 
finite differences for definition of set of operators 
|k = 2,A7 —1| at the second step is sufficient for 
stable calculations [2], 
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We remark, that no higher derivatives of Jacobian but only 
its degrees are used in suggested approach. It leads to 
significant simplification of code in the case of vector 
equation (1). So our approach is free of the main lack of 
existing Lax-Wendroff-type time integration schemes. 


IV. Numerical Results 


We construct FDSs with M={ 3, 5, 7, 9, 11} 
temporal and spatial order of accuracy for the numerical 
calculations. We use centered finite differences at the all 
steps ofFDS. We add the dissipative term of (M+l) spatial 
order of accuracy at the second step (11) for suppressing 
high-frequence oscillations in solution: 

v,=A"^W? (/,")) 02) 

We use well-known minmod-limiter\cite [6] for definition 
the value of a ( - + i/ 2 • ^ provides non-increasing of total 
variation of numerical solution. 

We checked stability of FDSs, used in this session. 
Spectral analysis in linear case and analysis of modified 
equation showed conditional stability of all used FDSs with 
classical Courant-Friedrichs-Lewy condition 

Cu = < i 

h 

Here 2 max is maximum velocity of disturbance 


propagation in each solved problem. We conduct all 
calculation with Cu = 0.5 in this session. 

We consider the accuracy test for nonlinear scalar Burgers 
equation [2]: 


8t 



= 0 . 


with initial condition /(0,x) = 0.5 + sin(^x) and 


2-periodical boundary condition. We conduct calculation up 
to time t ~ 0.5/ k , when solution is still smooth. Relations of 
L\ norm of approximation error to the number of grid point 
are presented in Fig. 1. 



Fig. 1. Approximation error of solution of Burgers equation 


Hereinafter k is estimation of order of approximation, 
defined by the calculation results. We can see that all FDSs 
aspire to their order of accuracy. We notice, that introduction 
of a dissipative term (12) reduces insignificantly an accuracy 
order. 


V. Concluding Remark 

We have developed a way to construct FDSs for PDE (1) 
with Lax-Wendroff-type time integration with (any) high 
order of temporal accuracy. This approach is based on 
classical Lax-Wendroff-type time integration scheme and one 
property of hyperbolic PDE (3), proved in this paper. Lor the 
first time LDSs with M ={ 3, 5, 7, 9, 11} temporal and 
spatial order of accuracy were constructed on the base of 
suggested approach. Numerical tests were conducted 
successfully to verify the accuracy for the case of the scalar 
hyperbolic equation. 
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